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I. INTRODUCTION 



The linear conductance of the a small metallic tunnel junction embedded in an electromagnetic 
environment of arbitrary impedance is determined in the semiclassical limit. Electron tunneling is 
treated beyond the orthodox theory of Coulomb blockade phenomena by means of a nonperturbative 
path integral approach. The frequency dependent conductance is obtained from Kubo's formula. 
The theoretical predictions are valid for high temperatures and/or for large tunneling conductance 
and are found to explain recent experimental data. 
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In recent years a great deal of experimental and theoretical work MM has explored Coulomb charging effects in 

systems with tunnel junctions. For small tunneling conductance Gt <C Gk, where Gk — e 2 /h is the conductance 

quantum, the effects are theoretically well understood in terms of the "orthodox" perturbative approach in the 

tunneling Hamiltonian |3|] . Present lithographic techniques make the region of moderate to large tunneling conductance 

experimentally accessible j|-|| . In fact when using Coulomb blockade devices in thermometry (7| or as highly sensitive 

electrometers [|3|, a large tunneling conductance is often desirable in view of the larger measuring signal. Several 

£j ' theoretical advances to describe strong tunneling have been made recently. These more sophisticated theories roughly 

split into two groups. One may use higher order perturbation theory in the tunneling Hamiltonian [|9|-[ll| and, based 

on this, a perturbative renormalization group approach |Fj, while the other set of papers starts from a formally exact 

path integral expressio n ||l3fl serving as a basis for analytical calculations |ld[-|l7|| , Monte Carlo simulations [|l8],[19| , or 

a variational approach j20|. Despite this progress a complete theoretical understanding is still missing. 

In this paper we focus on the frequency dependent conductance of a single metallic tunnel junction coupled to an 

electromagnetic environment of arbitrary impedance, and use a path integral formulation to analytically calculate the 

leading order quantum corrections. The semiclassical evaluation of the path integral is justified at high temperatures 

p^ ' and/or for large tunneling conductance and yields the frequency dependent linear conductance. So far experimental 

Q\ ', data p] are only available in the zero frequency limit where good agreement is found. 

The conductance of a tunnel junction is related to current fluctuations depending on the whole circuit. Via network 

transformations it is always possible to transform the linear electromagnetic environment into an admittance Y(lo) in 

series with a tunnel junction with tunneling conductance Gt and capacitance C biased by a voltage source V, cf. Fig 
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FIG. 1. Circuit diagram of a tunnel junction in series with an adnrittance. 

II. MODEL AND LINEAR CONDUCTANCE 

Our interest is in the zero bias differential conductance G = dI/dV\v=o where / is the current in the leads and V 
the applied voltage. Employing the Kubo formula 



G(w) = -L lim / dre^ (I(r)I(O)), 



(i) 



where the ^„ = 2im/Ti(3 are Matsubara frequencies, the conductance is related to the imaginary time current auto- 
correlation function which is most conveniently calculated as a variational derivative of a generating functional Z[£\ 
depending on an auxiliary field £(t). Introducing a phase variable tp, which is canonically conjugate to the charge 
transferred through the junction, the generating functional can be written as a path integral 
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with the effective Euclidean action 



Here 



describes Coulomb charging and 



S[<p,Z] = Sc[<p] + St[<p] + Sy[<p,Z] 
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quasi-particle tunneling across the junction. The kernel a(r) is determined by the tunneling conductance Gt and 
may be written as 
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where the Fourier coefficients are given by 






The last term in the action (0) describes the effective environment and includes the auxiliary field £(r) 
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where the kernel fc(r) can also be written as a Fourier series with coefficients 

ft?(KI), . 
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Here Y(s) denotes the Laplace transform of the environmental response function, cf. Ref. pip . Due to causality, for 
Re(s) > 0, one may write Y(s) — Y(is) where Y(u>) is the frequency dependent admittance of the environment. 
We now perform the functional derivatives explicitly and get for the correlator m& 



(I(r)I(O)) 
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where Z = Z[0] denotes the partition function and the current functional IUp, t] is given by 



I[<P,t] = *J dr'k(T-T'Mr'). 



(8) 



(9) 



The first term in (pfl is independent of the phase variable ip and can be handled exactly, whereas the second term 
cannot be evaluated without further approximations. Thus, we split the conductance in two parts G — G\ + G2. The 
first part stemming from the second order functional derivative of the action reads 



GiH = -^-^-k{-iu + 5)= Y(w), 
iTwj h 



(10) 



where the unique continuation |22| of k(oo) is obtained by defining the analytic continuation of the absolute value as 

z Re(z) > 



— z Im(z) < 



(11) 



The second part Gi of the conductance results from a product of two first order derivatives of the action, and we 
write 



G 2 (w) = -J-GjHu + tf) 

iriuj 



(12) 



with 



CaK) = | y £>M cxp {-J.%, 0]| Ffo>, i/„], (13) 

where the explicit form of the current functional (0) written in terms of the Fourier transform of the phase yields 
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(14) 



This is a formally exact representation of the linear conductance. 



III. SEMICLASSICAL LIMIT 



To proceed we expand the effective action around the classical path. Following the lines of Ref. |1G] we change to 
Fourier space diagonalizing the second order variational action. The eigenvalues are given by 



where 



h 2 p 

K v n) = —z-\v n \ G (y n )+Y(\v n \) 



G (v n ) = \u n \C + G 2 



(15) 



(16) 



describes the junction as a capacitance in parallel with an Ohmic resistor characterized by the classical tunneling 
resistance. Including the fourth order variational derivative of the action we get from (|13| ) the expression 
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(17) 



The convergence of this expansion depends crucially on the eigenvalues dTq). To estimate the validity, we consider 
the smallest eigenvalue in more appropriate units 
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(18) 



where j3 denotes the inverse temperature and Ec = e 2 /2C the charging energy. This eigenvalue has to be large 
compared to 1, and we see that the expansion is useful for large conductance Gt + Y(v\) 3> Gk and/or high 
temperatures /3Eq *C 2-7T 2 . We now perform the limit iv n — > w + iS. Consider the analytically continued eigenvalue 



h 2 p 



A(-tw) = -iu^t [G (u) + Y(lo)} , 



(19) 



where 



Gq(uj) = Gt — iwC 



(20) 



is the continuation of the Laplace transform Gq(—ilu + 5). The relative minus sign of the capacitive term is due to 
the usual definition of the Fourier transform in quantum mechanics, the electro-technical convention is obtained by 
replacing ui — > — ui. For small frequencies the continued eigenvalue (|l9|) is no longer large compared to 1 and we are 
faced with a problem of order reduction. This is handled systematically in Ref. fig] showing that inclusion of the 
fourth order variation is sufficient up to first order in 0Ec and Gk/{Gt + Y(fi)), respectively. In the reminder we 
omit the order symbol, but the meaning of the equations is always meant in a limiting sense up to this order. After 
performing the continuation we get for the conductance 



G{w) 
with an effective linear conductance of the junction 



G cS (m)Y{u) 



G eff (u>) =G T [1- U(uj)] - iuiC. 



(21) 



(22) 



This describes a linear element G*(lu) — Gt[1 — I4(u>)], depending on the whole circuit, in parallel with the junction 
capacitance C as depicted in Fig. 0a. Here 



un = _ 5> 
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(23) 



describes the suppression of the conductance due to discrete charge transfer. The general form ( |2lD is valid only to 
first order in j3Ec and Gk/{Gt + Y(vi)), respectively. A systematic treatment of higher order contributions does not 
allow for a description of the tunnel junction in terms of an effective linear element depending on the whole circuit. 
However, a partial resummation of higher order terms performed by a self consistent harmonic approximation EG] 
leads again to the form (pll) . 
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FIG. 2. Effective circuit diagrams for a tunnel junction in the semiclassical limit a) for arbitrary frequency and b) in the low 
frequency limit. 



IV. DISCUSSION 



For further discussion of the conductance and a comparison with experimental data we now restrict ourselves to 
ohmic dissipation Y(co) = Y. We then get for the effective linear element 



Gt 



= 1 - 



ip(l + u + iu) -V>(l+w) 4'{\ + u + uj) - ip(l + u)-\ j3E c 
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where i\) is the logarithmic derivative of the gamma function and 
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(24) 



(25) 



are auxiliary quantities. We also have introduced the dimensionless parallel conductance g = (Gt + Y)/Gk- The 
quantum corrections depend only on this combination of conductances. The real and imaginary parts of G*(uj)/Gt 
are depicted in Fig. for ftEc — 1 and various values of g. 
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FIG. 3. Real and imaginary parts of G*(u))/Gt in the ohmic damping case for (3Ec 
g — 1, 10, 40 and 100 in dependence on the dimensionless frequency Q — hu>/2irEc- 



1 and dimensionless conductance 



Due to the logarithmic behavior of the psi function for large arguments, the quantum corrections disappear non- 
analytically for large u. For small frequencies the effective element behaves like a renormalized admittance with an 
additional capacitance in parallel and we may define a renormalized capacitance C* and a renormalized conductance 
G* = G*(lu = 0), cf. Fig. pp. For the renormalized conductance we get 



G* 

G~, 



= 1 



which coincides with our previous result 
mental resistance. 
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L6[ . G* exhibits a nonanalytic behavior in the limit of vanishing environ- 
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FIG. 4. The linear conductance versus the dimensionless temperature for two dimensionless parallel conductance g — 4.2 
and 23.85 compared with experimental data (symbols) by Joyez et al. M. 

We compare our prediction (|2q ) with recent experimental data by Joyez et al. |J for dimensionless conductance 
g = 4.2 and 23.85. Fig. [| shows that in the limit of large conductance we are able to cover the whole temperature range 
explored experimentally with no adjustable parameter, whereas for moderate conductance only the high temperature 
part is covered by the semiclassical theory. 

The renormalized capacitance C* incorporates the linear part in u> of the imaginary part of G*(uj) and the geomet- 
rical capacitance C, and we get 
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(27) 



The correction shows a quadratic dependence on (3Ec- The renormalization is suppressed at high temperatures and 
also vanishes linearly for large conductance. 

In summary, we have derived an analytical expression for the frequency dependent linear conductance of a tunnel 
junction in the semiclassical limit. We have shown that this limit covers not only high temperatures but also large 
conductance and agrees with experimental findings. 
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